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Abstract 

We present a simple classical (random) signal model reproducing 
Born's rule. The crucial point of our approach is that the presence of 
detector's threshold and calibration procedure have to be treated not 
as simply experimental technicalities, but as the basic counterparts 
of the theoretical model. We call this approach threshold signal de- 
tection model (TSD). The experiment on coincidence detection which 
was done by Grangier in 1986 [20] played a crucial role in rejection 
of (semi-) classical field models in favor of quantum mechanics (QM): 
impossibility to resolve the wave-particle duality in favor of a purely 
wave model. QM predicts that the relative probability of coincidence 
detection, the coefficient g( 2 \0), is zero (for one photon states), but 
in (semi-)classical models g^(0) > 1. In TSD the coefficient g^(0) 
decreases as where £ d > is the detection threshold. Hence, 

by increasing this threshold an experimenter can make the coefficient 
g( 2 '(0) essentially less than 1. The TSD-prediction can be tested ex- 
perimentally in new Grangier type experiments presenting a detailed 
monitoring of dependence of the coefficient g^ 2 \0) on the detection 
threshold. 



1 Introduction 

We study the old problem of a possibility to construct a classical field model 
reproducing probabilistic predictions of quantum mechanics. The common 
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opinion, see Bell [Tj, is that this is impossible to do, but see [2] [EJ for 
numerous attempts to proceed with classical wave (oscillatory) prequantum 
models. The main argument is that for composite quantum systems (e.g., a 
pair of entangled photons) the correlation predicted by QM (and confirmed 
by experiment [I5]-[T7]) cannot be reproduced by models with (local) hid- 
den variables. This is also the common opinion that "Bell's theorem" pQ is 
formulated in so general abstract framework that it rejects all (local) models 
with hidden variables, including variables of classical field (e.g., electromag- 
netic) type. Although nowadays Bell's inequality and entangled systems are 
the hot topics [18], [19], the problem of a possibility of the classical field 
description of a single quantum system is also of the large importance for 
quantum theory, including quantum information theory. One of the most 
important tests of a possibility to represent a photon simply as a pulse of 
classical electromagnetic field is the experiment on coincidence detection in 
two output channels of the polarization beam splitter (PBS) [20], [21]. QM 
predicts that if the source can be considered as one-photon source, then the 
probability of coincidence detection equals to zero: one photon cannot be 
split between two channels, it is either in one or another channel. Any pulse 
of the classical electromagnetic field is split between two output channels of 
PBS. Hence, the probability of coincidence detection is nonzero. Of course, 
the real experimental situation is more complicated: for any "one-photon 
source" the probability of emission of e.g. two photons is nonzero (although 
it can be made very small), there is also the contribution of noise, including 
so called duck counts (i.e., counts in detectors in the absence of the source). 
Therefore, instead of the absolute probability of coincidence detection Pi 2 , 
experimenters use the relative probability [20j, [2"T] : 

S (2) (0) = -gk, (1) 



PiP 
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where Pu is the probability of coincidence detection in channels i — 1, 2 and 
Pi are probabilities of detection in corresponding channels. Known (semi-) 
classical models of QM predict that g^(0) > 1. If photon is not a classical 
pulse, one can expect that g^(0) < 1 (even by taking into account noise and 
emission of double photons). The first experiment of this type was performed 
by [2D] , see also [21] and [22] for review. It was shown that the number of 
double clicks is relatively small. This experiment played an important role 
in quantum foundations. As is commonly accepted, this experiment justified 
rejection of (semi-) classical field theories as "prequantum theories"; in par- 
ticular, photons cannot be interpreted as pulses of classical electromagnetic 
field. 

The aim of this paper is to show that the main reason of dis-matching of 
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predictions of (semi-) classical field models with QM and experimental data is 
that such models do not take into account the impact of detectors to creation 
of the quantum statistics. 

We present a classical field model which in combination with the pro- 
cedure of detection of random signals by threshold type detectors (which are 
properly calibrated) reproduces quantum probabilities and, in particular, the 
coefficient g^ 2 \0) << 1 for a sufficiently high detection threshold Hence, 
the prediction of our model differs crucially from the prediction of the known 
classical field models which do not take into account the evident fact that 
measurement is not simply detection (of the monitoring-type) of continuous 
classical signals, but creation of discrete counts with the aid of threshold 
detectors. Detectors also have to be properly calibrated. 

Our model, threshold signal detection model (TSD), predicts that the 
quantum prediction, g^ 2 \0) << 1, matches better measurements for high 
value of the detection threshold. The basic prediction of TSD is that the 
coefficient 

c d 

where K > is a constant depending on the signal. This K depends on the 
brightness of the source: higher brightness implies larger K. It also depends 
on elements of the density matrix p = (pij) corresponding to the prequantum 
signal. (In TSD p corresponds to the normalization (by the trace) of the 
covariance operator B of the prequantum random signal, i.e., p = B/TrB.) 
Larger p\2 implies larger K. 

It is interesting to compare this prediction with the real experiment. Un- 
fortunately, it seems that Grangier's type experiments with detailed monitor- 
ing of dependence of the coincidence probability on the value of the threshold 
have never been done; more specifically, that is: In Grangier's experiment, 
did the calibration of detectors play a crucial role to eliminate coincidences? 

This is clearly a crucial question. It is hard to give a definite answer 
as to how exactly it was influencing the result at a fundamental level be- 
cause Grangier et al [20], [2TJ didn't study the influence of threshold on the 
<7^(0) parameter he was measuring. It is nevertheless clear that with a lower 
threshold, he would have gotten a less good result. Let me translate the part 
of Grangier's thesis (21] where he explains for the first time the role of the 
threshold and how its level was chosen: 

"[...] In this configuration, the threshold has a double role of acquisition 
of timing information and of selection of the pulses (the too weak pulses 
are not taken into account). The problems connected to the choice of the 
discriminator threshold and of the high voltage of the photomultipliers are 
discussed in detail in reference [15]. We have in the present experiment 
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chosen a rather high threshold, which amount to give the priority of the 
stability of the counting rates and the reproducibility of the results, rather 
than to the global detection efficiencies." (I stressed with bold the important 
fact that Grangier proceeded with rather high threshold.) 

TSD provides a strong motivation to perform Grangier 's type experiment 
with monitoring of dependence of the coincidence probability on the detection 
threshold and source's brightness. 

TSD can be considered as measurement theory for recently developed 
prequantum classical statistical field theory, PCSFT, [23j,[24j. The latter re- 
produced all quantum averages and correlations including correlations for 
entangled quantum states. In particular, PCSFT correlations violated Bell's 
inequality. The main problem for matching of PCSFT and conventional QM 
was that PCSFT (nor other classical field models) was not able to describe 
probabilities of discrete clicks of detectors. In particular, PCSFT is theory 
of correlations of continuous signals. "Prequantum observables" are given 
by quadratic forms of signals. These forms are unbounded and this is not 
surprising that correlations of such observables can violate Bell's type in- 
equalities, see [22] for discussion and an elementary example. The condition 
of coincidence of ranges of values of quantum observables and corresponding 
"prequantum variables" plays a crucial role in Bell's argument. TSD solved 
the measurement problem of PCSFT. In the same way as in Bell's consid- 
eration, TSD operates with discrete observables. In particular, in the case 
of photon polarization (its projection to a fixed axis) TSD operates with 
dichotomous variables taking values ±1. 

TSD/PCSFT for composite quantum systems was presented in [26]. How- 
ever, random signals considered in [26J have a complex structure of temporal 
correlations, even in the case of a single system. In the present paper, we use 
simply a combination of the Wiener process (as the career of temporal cor- 
relations) and a Gaussian stationary process valued in the space of internal 
degrees of freedom of a quantum system (e.g., its polarization). We do not 
consider spatial degrees of freedom. 

To escape misunderstanding, we stress that the presented detection model 
is purely classical; in particular, detectors are classical detectors which are 
sensitive to the energy density of a signal in the domain of detection V. Such 
a detector clicks at the first instance of time r = t& when the total energy of 
signal in V approaches the detection threshold Ed- Since we consider random 
signals, this instant of detection is a random variable r = t(u>). Our main 
aim is to find its average f. The quantity 1/f determines the probability of 
detection. At the moment we proceed with such an operational description 
of the classical detectors. We plan to consider a more detailed scheme of the 
classical threshold detection in another paper in which we shall study classical 



4 



signals with spatial degrees of freedom, i.e., signals in physical space-time. 
(In the present paper we restrict the model to "internal degrees of freedom" 
such as polarization. The model with spatial degrees of freedom is essentially 
more complicated, since it involves processes with infinite-dimensional state 
space.) 

2 Threshold detection 

2.1 The class of random signals 

We consider a special class of classical random signals. The model is phe- 
menological: we cannot present physical motivations for selection of this 
class of signals, besides the fact that detection of such signals with the aid 
of threshold type detectors reproduces the correct quantum probabilities. 
Another class of classical random signals serving for the same aim was in- 
troduced in [26]. Signals considered in the present paper have essentially 
simpler temporal structure, simply the Wiener process. 

Structurally our model has some similarity with the prequantum model 
of Grossing et al [32]. Subquantum stochasticity is combined of the two 
counterparts: a stationary process in the space of internal degrees of freedom 
and the randmom walk type motion describing the temporal dynamics. 

We start with stationary signals. We proceed in the finite dimensional 
state space corresponding to internal degrees of freedom such as polarization. 
Generalization to spatial degrees of freedom is evident, but it has essentially 
more complicated mathematical structure. 

Let H be the m dimensional complex Hilbert space. Let = 4>{oj) be 
the if -valued Gaussian random variable with zero average and the covari- 
ance operator B. This operator is Hermitian and positively defined. Take 
in H an orthonormal basis {ej} and consider corresponding signal's compo- 
nents <f)j{uj) = (<p{u)), Cj). The mathematical expansion of the random vector 
4>{uj) with respect to the basis <f>(to) = J2j ( Pj( UJ ) e j physically corresponds to 
splitting of the signal into disjoint channels. We remark that correlations of 
signal's components are given by 

E(pi(j)j = (Bei,ej) = 6^. 

Now we introduce the temporal stochastics by simply using the one dimen- 
sional Wiener process w(t) which is independent from the stationary pro- 
cess 4>{oj). Consider the random (non-stationary) signal <f>(t,u) = w{t)(j){uj) 
and its components corresponding to internal degrees of freedom: <pj(t,u) = 
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w(t)(j)j(u). We remark that correlations of signal's components are given by 



min(t, s)bij. 



The energy of the ith component of the complex signal <fi(t) (at the instance 
of time t) is given by the square of its absolute value 



2.2 The scheme of threshold measurement 

We consider the measurement scheme in which each channel, i = 1,2, ...m, 
goes to a threshold type detector. We assume that all detectors have the same 
threshold > 0. The detection procedure under consideration is reduced to 
the condition of the energy level approaching the detection threshold. The 
instant of time r corresponding to the signal's detection ("click") by jth 
detector is determined by the condition: 



We remark that the instant of the signal detection is a random variable: 



Mathematically our aim is to find average of the instance of detection, f = 
Et. The quantity 1/f will be used to find the probability of detection, "how 
often the detector produces clicks," see section [3J 

We apply the mathematical expectation (average) operator to both sides 
of the detection condition (j2J) and we obtain 



Ei(t,u) = \&(t,cu)\ 2 . 



The total energy of the signal is 



E^u) = Y J \utM\ 2 = \m,cu)\\ 2 . 



Ej(r,u) = E d . 



(2) 



r = t(uj). 



EEj(r(cj),uj) = E d , 



(3) 



or 



Ew 2 (t(u),co)E\(P j (co)\ 2 = £ d . 
To find the first average, we use the formula of total probability: 



(4) 
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We know that, for the fixed r, Ew 2 (r,u) = r. Hence, 



Ew 2 (t(u),u) = / rfrP(r(w) = r)r = Et = t. 
Jo 

Thus, the detection condition fll]) has the form: 

fE\^{u)\ 2 = S d , (5) 

or 

rbu = £d- (6) 
We remark that r = Tj, i = 1, 2, m. Thus 

1 = ^- (7) 

3 Probabilities of clicks in detection channels 

Hence, during a long period of time T such a detector clicks iV click -times, 
where 

T bjjT , . 

iV,^ - = (8) 

To find the probability of detection and match the real detection scheme 
which is used in quantum experiments we have to use a proper normalization 
of Ni, This is an important point of our considerations. (The normalization 
problem is typically ignored in standard books on quantum foundations, 
cf., however, [23].) In QM-experiments probabilities are obtained through 
normalization corresponding to the sum of clicks in all detectors involved in 
the experiment, e.g., spin up and spin down detectors. 
Hence, the total number of clicks: 

We remark that the total number of clicks does not depend on the split of 
the signal into disjoint channels, i.e., on the selection of an orthonormal basis 
{e^} in H. In fact, 

£||0M|| 2 = TrB. 
The probability of detection for the jth detector is given by 

P t = Ni/N = (10) 
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In fact, this is the Born's rule of QM. Consider the operator 

p = B/TiB. (11) 

This is the Hermitian positive trace one operator; so it has all properties of 
the density operator used in QM to describe the state of a quantum system. 

Set Ci = \ei)(ei\, the orthogonal projection onto the vector e«. Then the 
equality for the probability of detection (JTU1) can be written as 

Pi = TrpQ. (12) 
This is the QM-rule for calculation of probabilities of detection. 



4 Coincidence detection 

Coincidence of clicks corresponds to matching of two conditions of threshold 
approaching corresponds to two constraints: 

£i(ti(u;),u;) = £ d , £ 2 (r 2 (u), u) = £ d = £d, (13) 

where matching has the form 

n(w) =r 2 (w) =r{u). (14) 

Our aim is to estimate the probability of coincidence P\2- To shorter notatin, 
we set £i{oS) = £j(r(u;); to) or even simply T i: i = 1,2. 

We consider the set of random parameters corresponding to coincidence 
detection: A 12 = {u : £j = £^i = 1,2}. (A 12 is the event of coincidence 
detection). We have to estimate its probability, P\ 2 = P{A\ 2 )). We shall get 
a rather rough estimate which, nevertheless, will be sufficient for our purpose. 
However, we shall see that better estimates of this probability will clarify 
essentially inter-relation between our "prequantum classical field theory", 
QM, and experiment. In principle, one may hope to derive an approximative 
expression for P\ 2 as we did for probabilities Pj,j = 1,2. However, this is a 
complicated probabilistic problem. 

We remark that A\ 2 is a subset of the set 

A lx2 = {u : £ ± £ 2 = £l}. 
Hence, P{A 12 ) < P(A ix2 ). And the set A x x2 is a subset of the set 

Ai x2 >q = {uj ■■ £\£ 2 > £d} 
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and hence P(A 12 ) < P(A lx2>£ 2). The latter probability we can (roughly) 
estimate by using Chebyshov inequality (which usage is standard for such 
estimates, cf. [27]). In the simplest form, for a random variable u = u(cu) 
and a constant k > 0, this inequality has the form: P(u : u > k) < 

'2 

ci- 



lia our case u = £i£ 2 ,k = £]. We have 



E£ x £ 2 



P(A 12 ) < (15) 



f 2 



We find this average by using the formula of total probability [?] 

POO 

E£ 1 {uj)e 2 {u) = / ES 1 {t,u)S 2 {t,u)P{t{u) = r)dr. (16) 



Thus our main problem is to find the correlation of two energies for each 
instant of time r. We have 

E£ 1 (t,u)£ 2 (t,u) = Ew\t)E\Mu)\ 2 \Mu)\ 2 - 

The first factor is known, Ew A (r) = 3r 2 . 

To find the second factor, we shall use general theory of Gaussian integrals 
on complex Hilbert space [28], see appendix. Consider in H (m-dimensional 
complex space) projection operators A k = | e^) (e^ | , A; = l,...,m. Set fA k ((p) = 
(Ak4>,(j)),cj) e H, the quadratic form corresponding to the operator A^. By 
fl27|) . appendix, we obtain 



Ef A J A2 = TrBA 1 TrBA 2 + TrBA 2 BA 1 , (17) 

where B is the covariance operator. We remark that 

by = PyTrB, (18) 

where p = (p^) is a density operator. So, we consider the prequantum ran- 
dom signal corresponding to the quantum state p and we study the problem 
of detection coincidence for such a signal 4>(t) = 4> p (t,cu). We have 

TrBAi = b ii: i= 1,2. 

In the same way 

TrBA 2 BAi = |( ei |5|e 2 )| 2 = |6 12 | 2 . 
Finally, we obtain (for the fixed instant of time r) 

ESfaufafau) = 3r 2 {b u b 22 + |6 12 | 2 ). 
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The formula of total probability ( 1161) implies 

E£ 1 {t(u);u)£ 2 (t{u)',u) = 3^(6 U 6 22 + |6 12 | 2 ), 

where r 2 = Et 2 . By ffTTl) we can rewrite this answer in terms of quantum 
mechanical density matrix 

££ x (t(u;);u;)£ 2 (t(u;);u;) = 3^TrB(p iiPjj + | Pij | 2 ), 

Thus by the Chebyshov inequality 

P(A 2 )<%^(puPn + \Pv\ 2 )- (19) 

The quantity 

A = (20) 

can be interpreted as the time parameter scaling the time intervals between 
coincidence clicks. Therefore £a = ATrB is average of signal's energy dis- 
tributed between clicks. Supppose that 

e = ^«l. (21) 

In this case 

P(A l2 ) < 3e 2 (p llPjj + \ Pij \ 2 ). (22) 

As usual [], set g^(0) = -j^p-, where Pi and P 2 are probabilities of detection 
in channels i and j, respectively (it is assumed that the later probabilities are 
positive). We proved, see sectional that in TSD model, these probabilities 
are equal to quantum probabilities: Pi = p n and P 2 = p 22 . Hence, 

S (2) (0) = ^<3e 2 (l + ^). (23) 

Since e ~ ^, by increasing the parameter £ d the coefficient #i 2 (0) can be 
done less than 1. Hence, our prequantum classical field model can violate 
(under the selection of a proper threshold) the inequality g^ 2 \0) > 1 which 
is valid for "standard" classical signal theory - its standard version which 
does not take into account the evident fact that measurement is not simply 
detection (of the monitoring-type) of continuous classical signals, but creation 
of discrete counts with the aid of threshold detectors. Detectors also have to 
be properly calibrated. 
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5 Appendix: Gaussian integrals 



Details of theory of integration with respect to Gaussian measures on complex 
Hilbert spaces can be found in [28J. 

Let W be a real Hilbert space. Consider a cr-additive Gaussian measure 
p on the a-field of Borel subsets of W. This measure is determined by its 
covariance operator B : W — >■ W and mean value m £ W. For example, B 
and m determine the Fourier transform of p : 

p( y ) = f e i{ *>*)dp(<j>) = e ^ By ' y)+i{m ' y) ,y £ W. 
Jw 

(In probability theory it is called the characteristic functional of the prob- 
ability distribution p.) In what follows we restrict our considerations to 
Gaussian measures with zero mean value: (m,y) = J w (y,ip)dp(ilj) = for 
any y £ W. Sometimes there will be used the symbol ps to denote the 
Gaussian measure with the covariance operator B and m — 0. We recall 
that the covariance operator B is defined by its bilinear form (Byi,y 2 ) = 

J(yiA)(y2A)dp((t)),yi,y2 £ W 

Let Q and P be two copies of a real Hilbert space. Let us consider their 
Cartesian product H = Q x P, "phase space," endowed with the symplectic 

operator J = ^ ^ ^ y ' ^ ons ^ er ^ ne c l ass °f Gaussian measures (with 

zero mean value) which are invariant with respect to the action of the oper- 
ator J; denote this class S(H). It is easy to show that p £ S(H) if and only 
if its covariance operator commutes with the symplectic operator, [23] . 

As always, we consider complexification of H (which will be denoted by 
the same symbol), H = Q © iP. The complex scalar product is denoted by 
the symbol (■,■). The space of bounded Hermitian operators acting in H is 
denoted by the symbol C S (H). 

We introduce the complex covariance operator of a measure p on the 
complex Hilbert space H : (Dyi,y 2 ) = J H (yi,(f>)((f),y2)dp((t>). Let p be a 
measure on the Cartesian product Hi x H 2 of two complex Hilbert spaces. 
Then its covariance operator has the block structure 

(24) 

where Da : Hi — )■ Hi and Dij : Hj — > Hi. The operator is Hermitian. Hence 
D ii = Ai, and -D* 2 = D 2 \. 

Let if be a complex Hilbert space and let A £ C S (H). We consider its 
quadratic form (which will play an important role in our further considera- 
tions) (ft — > = (A(j),(j)). We make a trivial, but ideologically important 
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remark: '■ H — > R, is a "usual function" which is defined point wise. We 
use the equality, see, e.g., [23J: 

/ f A ((/>)dp D (((>) = Tr DA (25) 
Jh 

Let p be a Gaussian measure of the class S(Hi x H 2 ) with the (complex) 
covariance operator D and let operators Ai belong to the class £ s (Hi),i = 
1, 2. Then 

/ fAMfAMdpty) = TrDuAt TiD 22 A 2 + Tr D 12 A 2 D 21 A 1 (26) 

JH 1 xH 2 

This equality is a consequence of the following general result [23]: 
Let p G S(H) with the (complex) covariance operator D and let Aj G 
C S (H). Then 

/ fA 1 {4>)fA 2 {4>)dp{4>), = TiDAiTiDA 2 + TrDA 2 -D^i- (27) 
iff 

Summary on coincidence probability: Prequantum classical field 
model with threshold and properly calibrated detectors violates predictions 
of standard classical and semiclassical models, cf. [J, and gives the predic- 
tion compatible with known experimental data. More detailed experiments of 
Grangier's type with monitoring of the dependence ofg^'(0) on the detection 
threshold and source's brightness are on demand. 
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